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Using the effective potential method, we computed one-loop corrections to the mass matrix of 
neutral Higgs bosons of the Non-Holomorphic Supersymmetric Standard Model (NHSSM) with 
explicit CP violation, where the radiative corrections due to the quarks and squarks of the third 
i—{ ' generation were taken into account. 

We observed that the non-holomorphic trilinear couplings can compete with the holomorphic ones 
■ in CP violating issues for the mass and mixing of the neutral Higgs bosons. 
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I. INTRODUCTION 
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In the Minimal Supersymmetric Extension of the Standard Model (MSSM) super-potential and soft breaking terms 
are generally considered as holomorphic functions. While the holomorphicity of the superpotential is obligatory for 
the MSSM, more generalized versions of the model may also include R-parity violating terms and/or non-holomorphic 
structures in the soft breaking sector of the theory [l], Q ■ 

Contrary to the SM where a unique Higgs boson resides, supersymmetric models predict extra Higgs bosons via 
introducing different Higgs doublets. With a number of supersymmetric models and various extensions, experimental 
I verification of Higgs bosons became one of the main objectives of the current colliders such as Tevatron [3j and LHC 
Compared with the SM, the allowed mass range of the lightest Higgs in the MSSM is somewhat more constrained. 
Indeed, the predictions related to the mass of the lightest Higgs in the MSSM give an upper bound m/, ~ 130 GeV [E[ 
, which can be fairly relaxed by certain extensions of the theory. In this sense, the Non-Holomorphic Supersymmetric 
Standard Model (NHSSM) requires the presence of the additional soft breaking parameters that can shift the upper 
0^ ! bound to a certain extend [|| which may be required in further Higgs searches. 

In addition to give a relaxation to the upper bound of lightest Higgs mass, these additional non-holomorphic soft 
breaking terms ensure extra degrees of freedom where, for instance, CP violating terms of the MSSM may get into 
trouble. It is explicitly shown in Q that the amount of CP violation present in the SM is not adequate to explain 
the observed baryon asymmetry in the universe, whereas, supersymmetric models offer novel sources of CP violating 
terms and especially Higgs interactions can play a key role in mediating CP violation. This issue is deeply scanned 
in the MSSM [8|. But, this should also be probed for the extensions of the minimal model where additional sources 
of CP violating terms exist. In this respect, the NHSSM is an interesting model with extra sources of CP violating 
terms in the soft breaking part of its Lagrangian. However before doing this, precise predictions are required for the 
Higgs sector of the model which does not exist in the literature. 

Hence, in this work, our interest focused on the neutral Higgs sector of the Non-Holomorphic Supersymmetric 
Standard Model (NHSSM) with R-parity conservation and explicit CP violation. We assumed that CP is explicitly 
violated in the Higgs sector of the NHSSM and looked for its impact on the mass and mixing of the neutral Higgs 
bosons, which may be important for the near future. The possibility of non-holomorphic structures in the soft breaking 
is realized in the literature. For a detailed list of issues ranging from b — > 37 decay to the Renormalization Group 
Equations (RGEs) of the non-holomorphic supersymmetric model, we refer to [9(, a possible explanation for the source 
of the NH structures can be found in [l(| • 

The rest of the paper is as follows: In the following section we first described the basic low energy structure of the 
NHSSM. Analytical results for the mass matrix of the neutral Higgs bosons of the NHSSM are derived in the same 
section where the one- loop CP violating effective potential is calculated by considering only top and bottom sectors. 
Section III is devoted to numerical analysis where the impact of non-holomorphic trilinear terms are investigated. 
And we concluded in section IV. 
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II. NHSSM 

In general the supersymmetry (SUSY) breaking sector is parameterized via holomorphic operators which must be 
soft i.e. the quadratic divergences must not be regenerated [ll|. However, the MSSM can be extended via introducing 
new soft operators including R- violating and/or non-holomorphic terms in the soft breaking sector of the theory. In 
this sense there are different non-holomorphic models based on different approaches (see [l[ and 0]). Here we follow 
the easiest path in which R-parity violating terms are ignored and the problematic Higgsino mass term (pi) is absent in 
the superpotential. Under these assumptions, the NH version of the minimal supersymmetric model can be described 
by the superpotential 

W = Q-H u Y u U-Q-H d Y d D-L-H d Y e E (1) 

where our conventions are such that, for instance, Q -H u = Q T (ia^Hu = tijQ l W u with ei 2 = —£21 = 1- In the MSSM 
the breakdown of supersymmetry is parameterized by a number of holomorphic soft operators [Til ] 

- Csoft = Q f in|g + WmfjU + &m 2 b D + Dm\L + E^m%E + - (M 3 A?A? + M 2 A^A^ + M X X § \ § + h.c.) 

+ m 2 Hu HlH u + m\H\H d + (m 2 3 H u ■ H d + h.c.) + (q ■ H U Y*U - Q ■ H d Y%D - L ■ H d Y*E + h.c.) .(2) 

Here mj~ ^ are the soft mass-squareds of the scalar fermions, Y^ d e are their associated holomorphic trilinear 

couplings, and finally, M\, M2, M3 are, respectively, the masses of hypercharge, isospin and color gauginos. For the 
description of the Higgs sector soft masses m 2 H , m 2 Hd and m 2 are used. 

In the MSSM one can introduce the CP violation through the Higgs superpotential and the soft supersymmetry 
breaking terms, however, as has been shown explicitly in [l2.ll3j]. in supersymmetric theories which do not have pure 
gauge singlets in their particle spectrum, the holomorphic supersymmetry breaking terms do not necessarily represent 
the most general set of soft-breaking operators. Indeed, for instance, the MSSM spectrum does not consist of any 
gauge singlet superfield, and thus, its soft breaking sector must necessarily include the following soft breaking terms 

C' soft = n'H u ■ H d + Q H^Y'tO + Q H^Y'^D + L H^Y'^E + h.c. (3) 

in addition to those in ([2]). Here Y'^ d e are non-holomorphic trilinear couplings which do not need to bear any 
relationship to the holomorphic ones Y^ d in ([2]). The only exception to this can be imagined as a unique common 
term at very high scales but even in that case, due to renormalization group running effects it is good to assume those 
new trilinear couplings to be completely different from the ordinary ones. Since these non-holomorphic couplings are 
perfectly soft they must be taken into account when confronting the MSSM predictions with experimental data. 

For possible variants of the non-holomorphic model notice that the original fj, term can be protected in the super- 
potential, in this case, the soft breaking // can stand alone or replaced with // — /1, then the m\ term of the soft 
sector can be written as m| = B(/j, — pi'). But we have chosen to deal with only one /i parameter for which the prime 
symbol will be dropped from now on. 

A. Analytical Results for the Neutral Higgs Bosons 

As in the MSSM, the CP violation mixes neutral Higgs bosons and hence they will be depicted as physical mass 
eigenstates ^1,^2 and /13. The classical potential for the neutral Higgs fields can be written for the NHSSM like in 
the MSSM [3 as follows 

V = m ^Hl\ 2 + m 2 Hd \H°\ 2 -{ m lHlH° (4) 

Here g 2 = g\ + gy and AV refers to loop corrections which will be computed using the effective potential formalism 
and the symbols 52 and gy stand for the SU(2) and U(l) gauge couplings respectively. Contrary to the MSSM, in (|4|) 
|^| 2 contributions coming from F— terms are absent which stem from the superpotential (TTJ). In (j4|) we allowed m 2 to 
be complex by assuming: 



m 2 = \m 2 \e v 



(5) 



3 



This phase $ can be set to zero at the tree level, but due to loop corrections it should be protected. Our results 
for Higgs masses will depend on this and another phase coming from one of the Higgs doublets. Now, the neutral 
components of the Higgs doublets can be expanded around their vacuum expectation values (vevs) as 

1 e ie 
H d = -j={v d + <f) d + ia 1 ) , H u = -= (v u + 4> u + ia 2 ) (6) 

in which v 2 = v 2 + v d — (246 GeV) 2 . Using the ratio of vevs we define tan/3 = v u /v d . In the above expression, a 
phase shift e is attached to neutral part of th e up Higgs doublet and this phase should be fixed by true vacuum 
conditions considering loop effects (see [HI, [HI, Il7l . for details) . 

It is important to emphasize that without loop corrections the tree level Higgs potential ^ predicts the lightest 
Higgs boson to be lighter than the Z boson. Hence, as in the MSSM, sizeable radiative corrections are also needed in 
the NHSSM to satisfy the LEP bound of m/, ~ 114 GeV. In this work we did not consider LEP excess events which 
indicate on a possibility of even a lighter Higgs [l9j . however, the excess events may still remain as another issue to 
be considered within the context of the NHSSM. 

In the MSSM, the radiative corrections [IHGil are dominated by loops of the top (s)quark, and to a lesser extent, 
by those of the bottom (s)quark, tau (s)lepton, charginos and neutralinos [17j . A particularly useful framework for 
computing the radiative corrections in the Higgs sector is effective potential approach [2(J. At the one- loop level we 
can write the contributions of all the relevant particles (coupled to Higgs bosons) as 



64 n z 



i / M. 3 
M (ln A^-2 



(7) 



Here, A is the renormalization scale and M. is the field-dependent mass matrix of quarks and squarks with overall 
factors -12 and 6, respectively. The additional contributions coming from charginos, neutralinos, etc. are ignored in 
this work. 

To proceed, the relevant fermion masses should be stated in a field dependent manner as in the background formed 
by the neutral components of the Higgs fields, for instance, the squared-mass of bottom and top quarks are given by 

ml = \h b \ 2 | H° | 2 , m 2 = N 2 | Hi | 2 (8) 

and those of the scalar quarks are 

M2 _ ( ml + ml + ^ + g 2 Y K\HX - \H° d \ 2 ) h (A b H° A' b * iff) \ 

b ~\ K ( A t H T ~ A 'b H u) rnl+ml + yU\H° u \ 2 -\H d \ 2 ) ' [9 > 



4 " { K (A$H° - A'H°) m? + ml - \gl{\H« u \ 2 - |ff«| 2 ) J ' {W) 



^ mt - ^(3 ff 2 - g 2 Y )(\Hl\ 2 - |H0| 2 ) h t (A t H° - A> t *Hf 



In writing the squark mass-squared matrices we have introduced some notationally simplifying definitions such as 
(mj~)33 = m 2 ^ and (Y'^ ) 33 = h b A\. In ([9]) and p0|) . the main effect of non-holomorphic trilinear couplings is 
to replace the /i parameter in the holomorphic MSSM in a flavor-dependent way and this shift alone tells us that 
the fi parameter seen by Higgsinos is completely different than what is felt by the scalar fermions. It is possible to 
back-transform A' t ,A' b — > /x to obtain the MSSM results, but the reverse is not true. In other words, the indirect 
relation between scalar fermions and charginos or neutralinos over the fi parameter is completely vanished in this 
version of the NHSSM. From now on, bounds on the /i parameter (i.e. obtained from charginos) have no restriction 
on scalar fermions anymore. Besides this the mass of the Higgsinos is the same with the MSSM if we assume /i as an 
input parameter. 

In this work, rather than providing a general analysis of Y' u d e in regard to MSSM phenomenology (see [2l| 
attempts in this direction) , we will focus mainly on their influence on Higgs- fermion-fermion couplings (especially for 
hi — > bb decay) in order to determine their distinctive features and observability in collider experiments. Concerning 
this class of observables, the primary objective would be to determine sensitivities of Higgs boson masses and mixings 
to the non-holomorphic couplings Y' u d . For this purpose, to leading order, it suffices to consider only the top and 
bottom quark sector. 

Now, for later convenience, we introduce E and A symbols such that (s)top and (s)bottom mass eigenvalues can 
be written simply as 

= fs B T VA^) /4 (11) 
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satisfying hierarchical order m 2 ~ 
to field dependent m% 

Jl,2 



< m% with f 

h 



b, t. Of course, the sfermion masses appearing in correspond 
evaluated in the electroweak vacuum. The explicit form of our definitions can be found in 

J ! 1 

the Appendix. 

The calculation proceeds by plugging the field-dependent eigenvalues into the potential. The mass matrix of the 
Higgs bosons is given by the second derivatives of the potential (at vanishing external momentum) . For this aim, the 
minimum of the potential should be obtained, which can be extracted from the first derivatives of the potential V (|4|). 
In turn, m 2 H ^ , m 2 Hd and 7713 can be expressed in terms of functions of the parameters appearing in the loop-corrected 
Higgs potential. We define 



T t = dV/dV, 



(12) 



for the first and second derivatives respectively with Vl^, ^fj = <f> u , 4>d, oi, 02 and both T and Ai 2 should be evaluated 
at vacuum conditions \&j = 0. Among the stationary relations 7^ u and T,p d are linearly independent. But T ai and 
T a2 can be expressed in terms of each other. Hence we can express m 2 as follows 



16-7T 
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2rn? ln-A 
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2mr In 



A 2 



(13) 



In this equation 2/ = (AfA'^e 10 ) describes the amount of CP violation in the sfermion mass matrices. Notice that the 
combination of phases 9 + $ is re-phasing invariant and validity of (fTB"]) can be checked from [12] in the (A' t , A' b ) — ► fi 
limit. During the numerical analysis we fixed 9 = —n/2 and determined $ in accordance with the input parameters. 

After obtaining true tadpoles correctly, Higgs boson mass-squared matrix (M 2 j) is acquired in the base of 
{(/i„, fa, 01, 02} in the form of a symmetrical 4x4 matrix. The eigenvalues of this symmetric mass-squared ma- 
trix correspond to a massless Goldstone boson and three physical neutral Higgses (m 2 li ,mf l2 and rn\^j. They can be 
used for numerical purposes, but for analytical purposes it is useful to perform the following unitary transformation; 



S T M 2 S where S 



1 
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and 77 



C{3 -Sp 



(14) 



This transformation allows us to redefine A4 2 as a symmetric 3x3 matrix in the basis {4> Ul 4>di a} where a is defined 
as a linear combination of a\ and a-i (a — sin/3ai + cos j3 02). For instance component of the redefined 

matrix becomes, 



M 2 



33 



VdV u 



32ir 2 VdV u 



In- 



,w 2 cos(6> + $) 3\h t \ 2 v 2 Kt, m \ m \ 



A 4 



3\hb\ 2 v 2 K blii m i m i 



32n 2 VdVu 



A 4 



3|/i t | V£ T (8\h t \ 2 v d v u l 2 - A T TZt 



32n 2 VdV u A 1 ( 



3/2 



■In 



2>\h b \ 2 v 2 Y, B (S\h b \ 2 v d v u ll - A B U b 



2 

mi 

t2 



32-K 2 VdV u A^ 



3/2 



■In 



2 

t>2 



■—^ — r- {l6\h b \ 4 v 4 v d v u A T l 2 - A B ((v 4 d + vi)A T (2\h b \ 2 K b + 2\h t \ 2 K t ) 

+Av d v u (-4\h t \ 4 v 4 l 2 +v d v u A T (\h b \ 2 K b + \h t \ 2 K t )))} 



(15) 



The explicit form of the symbols given here can be read from the Appendix. We refer to the same place for the 
remaining five entries of the symmetric Ai matrix. 



III. NUMERICAL ANALYSIS 



In this part, based on our analytical results, our aim is to show how the interplay of the non-holomorphic couplings 
with holomorphic ones can change the mass and the mixing of the neutral Higgs bosons. 

During the analysis, to respect the collider bounds, we require the scalar fermion masses satisfying mj > 
100 GcV and generally our results cover the LEP bound m^ ~ 114 GeV. Our basic input parameters are 
Mai mQ^mjjjmf,, A t ,A b ,A' t , A' b and tan/3. During the analysis we fixed 9 = — 7r/2,A = 0.5 TcV and determined $ in 

accordance with the input parameters. The true phase of CP violation can be defined as 9 e ff — arg{AfA'^e ) for top 
and bottom sectors and hence variation of $ is never presented. Instead we concentrated on the trilinear couplings 
and looked mainly for the mass and the mixings of the Higgs bosons under CP violating non-holomorphic trilinear 
couplings. 
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In the numerical analysis we have taken Ma as one of the input parameters and forced Ma = .M33, which is 
slightly different from the selection of ref. [12]. Alternatively, mass of the charged Higgs boson can be used as an 
input parameter as it is usually done in the MSSM literature (i.e. see 8}). To make the neutral Higgs mass matrix 
diagonal, we also defined an orthogonal matrix O such that 



diag{m 2 hi ,m 2 h2 ,m 2 ha ) = T M 2 Q, 



(16) 



then, an additional parameter can be defined to describe the CP composition of the neutral Higgs bosons as [22 



\0, 



|O l2 | 2 



\o. 



(17) 



With these definitions CP composition of neutral Higgses under the influence of non-holomorphic trilinear couplings 
can be presented. Additionally, we have selected the decay width of neutral Higgs bosons into bb as a testing ground. 
Of course (hi — > bb) has massive background, but according to the SM up to Mh ~ 130 GeV this channel is the 
dominant one. So we have chosen this channel for its simplicity to show the impact of the phases. Notice in the SM 
that the width of this decay is 0.0035 to 0.086 GeV for m h between 120 - 160 GeV ;23]. But the NHSSM has potential 
to change it sizably, as we will see. The partial decay width of a neutral Higgs boson hi into a pair of bottom and 
anti-bottom quarks is given as _23] 



r(hi-66) = 



i 9 W b m ht 
32Trm 2 v 



'1- 



o 2 



cos 2 



(l--^)+tan 2 /?0 2 3 



(18) 
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FIG. 1: The lightest neutral Higgs mass mux (left panel) and its CP violating parameter qi (right panel) of the MSSM (solid 
lines) and of the NHSSM (dashed lines), against m,Q for Ma=1 TeV (thick lines) and Ma=130 GeV (thin lines). Inputs: 
fj, — 500 GeV, A' t = 5m,Q, A' b is fixed at 0, trig = rtij} = tuq, A t — A b — 2ttiq and m,Q scans from 0.3 TeV to 1 TeV, tan/3 is 
fixed at 10. 



Now let us present our numerical results. We start with the difference of in the MSSM and in the NHSSM 
for CP violating case (see @ for CP conserving case). For the model under concern, the /i parameter of the NHSSM 
is responsible for higgsino masses as given in the soft breaking sector ©, and can be bounded from chargino and 
neutralino masses as in the MSSM. However in the NHSSM, // does not exist in the sfermions, so whenever we assume 
A'j 7^ /i we are considering the NHSSM. Thus for the differently selected fj, and A' t values model dependent effects can 
be observed on the mass and the mixing of the Higgs bosons. 

In this sense, Fig. [1] depicts the mentioned differences of the MSSM with fixed /i and of the NHSSM with different 
A' t and A' b values. The mass difference, as can be seen from the left panel of the figure, can be around few MeV or a 
few GeV and increases as the squark mass increases. This mass difference is determined by the magnitude of Ma and 
is more visible when Ma is close to m^. A similar observation can be extracted from the CP violating parameter of 
the lightest Higgs boson (right panel). Again when Ma ~ mh ± , thisparameter shows that the CP composition of the 
lightest Higgs can be enhanced as should be expected in the MSSM [8J|. In comparison to the MSSM, the CP violating 
parameter of the NHSSM can be smaller or larger than that of the MSSM predictions, which is again determined by 
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Ma- For instance, for Ma = 130 GeV the CP violating parameter takes a\ ~ 0.4 which is approximately half of 
the MSSM's prediction but for Ma — 1 TeV we observe af H ~ 0.25 which is approximately two times larger than 
the prediction of the MSSM. Notice that using this parameter one can determine whether hi will behave similar to 
the SM's Higgs boson or not. In sum, the effects of the extra parameter of the NHSSM (A' t ) is more visible when Ma 
is close to nih- This is something predictable because in this case large mixing occurs between the would-be CP-odd 
and CP-even Higgs bosons. On the other hand, when Ma is large it may be hard to effect the mass of the lightest 
Higgs boson with the NH terms. But even in this case its coupling could be very different from the prediction of the 
MSSM, which can be read from the right panel of the figure. 

In the same figure, the relaxation on the bound of sfermion masses can also be deduced for the NHSSM, i.e. as can 
be seen from the left panel of Fig. Q] for m hl > 115 GeV and Ma = 130 GeV the MSSM demands m squar k s > 450 
GeV, but in the NHSSM this bound relaxes to m squar k s > 350 GeV, for the selected range of parameters. 




arg(AJ) arg(Ai) arg(AQ 



FIG. 2: Masses of all neutral Higgs bosons rrihi 2 3 (l e h)> their CP-violating mixing angles 01,2,3 (center) and their decay widths 
into bb pair (right) versus the argument of the top trilinear coupling arg(A£). The dimensionful terms are given in GeV. Inputs: 
tan/3 = 10, rriQ = = rrif, = 1 TeV and At = At = \A' t \ = A' b = 2m,Q, Ma = 200 GeV. Note that for this and latter figures 
line formats are as follows, solid line (hi), dotted line (/12) and dashed line (/13). 

In Fig. [2] we present the phase dependencies of the masses of neutral Higgs bosons hi , hi and h% in left panel, 
CP violating parameters (ccj) for each of the mentioned bosons in middle panel and the corresponding decay widths 
T(hi — > bb) in right panel, against varying phase of the non-holomorphic trilinear coupling A' t . As can be seen from 
the first panel of Fig. [2] all of neutral Higgs bosons are sensitive to the phase of the non-holomorphic trilinear coupling 
A' t , with varying order. As a result of this phase the lightest Higgs boson can be made completely CP-odd or CP-even. 
One can easily recognize from the middle panel of Fig. [3 that CP violating mixing angles (ai) are suppressed for 
arg(AJ) ~ 7r/2 and have a sharp maximum value for arg(A(.) ~ 71-/3 and ~ 27r/3. Additionally the partial decay 
widths of the neutral Higgs bosons are very sensitive to the phase of A' t , can exceed the prediction of the SM for each 
of the bosons, i.e. T(hi -^bb)< 0.14 GeV is possible in the NHSSM as can be seen from the right panel of Fig. [2l 




arg(A^) arg(A(,) arg(A{,) 



FIG. 3: The same with Fig. [2] but now tan/3 = 50, Ma = 130 GeV, plots are presented against varying A' b . Inputs: 
A t = A b = A' t = \A' b \ = 2m Ci 
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In order to show the importance of A' b contribution, we present Fig. [3l in which all the input parameters are the 
same with Fig. [5J but now tan (3=50, Ma = 130 GeV. Here it is interesting to observe that the coupling of the lightest 
neutral Higgs boson is very strong, sensitive to the argument of the A' b , additionally T(hi — > bb) can be as large as 
>3 GeV (solid line of the right panel) which is well above the SM prediction. A common property of the Figs. I2l3l is 
that when 6 e ff ^ 0, mass difference of hi and /13 bosons increases. 

Notice that if were to assign arg(Af) = —arg(A^) then there would be no variation for the masses, the CP 
compositions and the partial decay widths of Higgs bosons which can be important for CP violating issues such as 
Electric Dipole Moments (EDMs) of fundamental fermions. This can be seen from the definition of the CP violating 
parameter If. 

It can be inferred from the presented figures that the NH soft breaking terms can yield sizable variations on the 
masses, CP compositions and partial decay widths of neutral Higgs bosons. These decay widths covers the range from 
~ SM values up to 3.5 GeV, CP violating parameters can be obtained from zero to one, thanks to the NH terms. 
Some of these results can also be simulated with a complex /x parameter of the MSSM (when A' t = A' b = fi), but 
the option of replacing this parameter with the non-holomorphic ones should be seen as an attractive alternative for 
the continuing and coming Higgs searches. While not presented here, it is easy to guess that the couplings of Higgs 
bosons to vector bosons are also sensitive to the mentioned NH terms. 



IV. CONCLUSION 

In this paper, we studied the mass matrix of the neutral Higgs bosons in the NHSSM with explicit CP violation at 
the one-loop level. For doing this, we first obtained analytical expressions for the mass matrix of neutral Higgs bosons 
(see the Appendix) and performed a numerical study based on the new sources of CP violating trilincar terms. 

In order to maximize the impact of the NH terms, among many possible parametrizations of the NH model, a special 
one is selected in which the indirect relation between scalar fermions and inos over the /i parameter is disappeared. In 
this version of the NH model, \i term is absent in the superpotential and it exists in the soft breaking part as a mass 
term for Higgsinos. We observed that this can heavily effect the mass and the mixings of the neutral Higgs bosons of 
the MSSM. 

During the numerical analysis we intentionally considered beyond the MSSM scenarios (such as A' t ^ A' b ) and 
observed that not only holomorphic soft breaking terms but also non-holomorphic terms can induce sizable amount of 
CP violation in the Higgs sector. In order to show this we studied numerically partial decay widths and CP-violating 
parameters of all the neutral Higgs bosons. We believe the selected ranges of our examples can be important for 
the continuing and upcoming (i.e. see [24]) Higgs searches with a generalized soft breaking MSSM. Additionally, this 
issue should be probed deeper because new sources of CP violating terms consisting with current collider bounds can 
be useful, for instance, to relax the electron and neutron EDM bounds on the CP violating terms of the MSSM [25| . 
Analysis of various observables ranging from b — * 57 decay to EDM constrains can shed further light on the structure 
of non-holomorphic models. 
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APPENDIX A: DEFINITIONS AND MATRIX ELEMENTS 

The definitions that appear in our calculations are as follows: For the proper treatment of the amount of CP violation 
we collected imaginary and real parts of frequently appearing terms as T b = Im(Af, A' b e 10 ), X t = lva(A t A' t e l9 ), similarly 
TZb = Ke(Ab A' b e l6 ) and TZt = Ke(A t A' t e l9 ). This enables one to decompose scalar fermions into parts as 

S B = 2m\ + 2m\ R + v 2 u X Gb + v\ (2\h h \ 2 ~ £ Gb ) , A B = n\/A 2 Gb + 8\h b \ 2 (\A b \ 2 v 2 d + v u (\A' b \ 2 v u - 2v d TZ b )) (Al) 
S T = 2m 2 L + 2mf fl - v 2 H Gt + v 2 u {2\h t \ 2 + £ Gt ) , A T = 4/A 2 Gt + 8\h t \ 2 (\A' t \ 2 v 2 + v u (\A t \ 2 v u - 2v d K t )) (A2) 
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where 

ki = A Gb (2ml R - 2m\ L + (v d 2 - v u 2 ) A Gb ) , k 2 = A Gt (-2m| L + 2mf R + (V - w u 2 ) A Gt ) (A3) 

k 3 = A Gb (2m? fl - 2m 2 L + (v d 2 - iv u 2 ) A Gb ) , k 4 = A Gt (-2m| L + 2m| R + (V - 3i;„ 2 ) A Gt ) (A4) 

k 5 - A Gb (2m 2 fl - 2m\ L + (3v d 2 - i> u 2 ) A Gb ) , n e = A Gt (-2m 2 l + 2m\ R + (3v d 2 - v u 2 ) A Gt ) (A5) 
A Gt = (-3. 9 2 + 5g 2 Y ) /12, S G( = - (.g 2 + .g 2 ,) /4, A Gb = (3.g 2 - .g 2 ,) /12, £ Gb = (.g 2 + g 2 -) /4 (A6) 

Additionally, we defined the following quantities such that entries of the Higgs matrix can be expressed in simpler 
forms. 

Xi = 4 ({4A b 2 \h b \ 2 + Kl ) v d - A\h b \ 2 v u TZ b ) , X 2 = 4 ((^A' t 2 \h t \ 2 + k 2 ) v d - 4\h t \ 2 v u n t ^ , X 3 = 4 (4A b 2 \h b \ 2 + k 6 ) 

X4 = 4 (4A;> t | 2 + k 6 ) , X 5 = 4 (4A' 6 2 |/i6| 2 ^ - «i«u - 4|/i fc | 2 i; ( i^ h ) , xe = 4 (4A t 2 |/i t | 2 w u - k 2 w u - 4|fc t | 2 i; d ^ t ) 

X7 = -8 (v d v u A 2 Gb + 2\h b \ 2 TZ b ) , X s = -8 (v d v u A 2 Gt + 2\h t \ 2 TZ t ) , X 9 = 16A' b 2 \h b \ 2 - 4k 3 , X io = l&A t 2 \h t \ 2 - 4k 4 . 

Using the definitions given above, the elements for the mass matrix of the neutral Higgs bosons due to radiative 
contribution of quarks and squarks are obtained as follows: 



9 2 v d 2 2 v u cos[e + *] 3 (8^ d 2 X2 A T S Gt + 2 X i A T S T + v d ( X2 2 - 2x 4 A T ) Sr) . m u 



4 v d 10247r 2 t; d A 3/2 

3\h b \W. ln m| 3 ( 2 (Xi + X 2 )A B A T + u d (x 2 2 A B + ( X i 2 - 2( X3 + X 4 )A B ) A T )) 



1-2 



8tt 2 A 2 5127r 2 v d A s A T 



3 (2 X iA B + v d (xi 2 - 2 X3 Ab)) S b + 24^ 2 X iAb {-2\h b \ 2 + S Gh ) ^ ^ 

1024^ 2 w d A 3 B /2 11 



*B 62 

3(^ 3 (A Gt +S Gt )+2NV,ft t ) m 2 m 2 2 'W (a 2 ,, + (-2|fe b | 2 + E Gb ) 2 ) + 6|fe b | 2 ^ b m 2 ^ 
647r 2 w d n A 4 + 647r 2 w d n A 4 



4/( 2 !-2 ,2/ M , *n i 3 (X2X6A B + (xiX5 - 2(x7 + X8)A B )A T ) 

•^12 = -4.9 v d v u + m 3 (- cos[6» + $]) + 5127r 2 A B A r ^ ^ 

3 (~8|fe t | 2 ^X2A T - 4^x2A T S Gt + 4v d x6A T Z Gt + X2Xg^t - 2x 8 A T S T ) ^ 

10247r 2 Ay /2 "if 2 
| 3 (-8|fe fc | 2 ^X5A B + XiXsgg - 2x 7 AbS b - 4^xiA B E Gb + 4v d x 5 A B Z Gb ) ^ 

1024tt 2 A| /2 11 ™ 2 2 

2 ™2 ^_ .. v, I „it 12 , v. \\ m 2 m 2 



| 3 ( X 8 ~ 8^ M S Gt (2|fe t | 2 + S Gt )) ^ mfmf a | 3 ( X7 - 8^ M S Gb (-2|fe b | 2 + S G „)) ^ 



61 6 2 



512tt 2 A 4 512tt 2 A 4 



2 _ 2 v d Sm[0 + $] 3(\h b \ 2 (v 2 xi-2v d A B )A T l b + \h t \ 2 A B (v 2 X 2-2v d A T )l t ) 
• Ml3 ~ ma S + 32^A^ (A9) 

j\h t \ 2 v d l t 3\h t \ 2 (-2v d A T X T + v 2 (4v d A T Z Gt +X2Z T ))lt m j 

32^v n A 4 + 64rr 2 «Af 



3|fe b | 2 ^J b 3|fe fc | 2 (-^ 2 XiSb + 2^A B S B - 4^Ag (-2|/i b | 2 + S Gb )) 1, 

32^v n A 4 647r 2 T,A 3 / 2 



9 

ln^, 

m~ 

B b 2 



9 



Ml 3 = m 2 



■j u Sin[0 + <E] 3 (|/t & | 2 (^ 2 X 5 ~ 2vuAb) A T l b + |/t t | 2 A B (v 2 X 6 - 2v u A T )l t ) 



3\h t \ 2 (4v 2 v u A T (2\h t \ 2 + E Gt ) - v 2 X e^T + 2u u A T S T ) J t | 3|/i t | 2 u u J t ^ m^m^ 



32 7 r 2 -t;A i 3A T 



64ir 2 vA 3 T /2 



2 

*2 



-2v u A b ^b + v 2 { X 5^b - Av u A B Y, Gb )) l b 



64n 2 vA 



3/2 



In- 



3\h b \ 2 v u l b ™* ™t 



32ir 2 v 



32ir 2 v 
■In 



■ In ■ 



A 4 



A 4 



(A10) 



KA i _9 2 v u 2 . 2 v d cos[6 + <S>] 3(2(x5+X6)A B A T + Utl (x6 2 A s + (x5 2 -2(xio+X9)A B ) A T )) 

AI22 — ^ r m? I — — — 5 (All I 

4 J v u 5127r 2 w M A B A T 

3|/i t | 4 « u 2 m 2 3(2x 5 A s S B +t; u (x5 2 -2x9A B )S i3 -8t; M 2 X5A B I] Gi) ) ^ 



8tt 2 A 2 10247r 2 « u A| /2 ™| 

-2Av u 2 X qA t (2\h t \ 2 + S Gt ) + 3 (2x 6 A T + ( X6 2 - 2 X iqA t )) S t mg i 

1024^ U A?/ 2 n ml 



3v u 3 (a^ t + (2N 2 + E Gf ) 2 ) + 6\h t \ 2 v d K t m » m a 3 ( Wu 3 ( A 2 ^ + S 2 J + 21^12^. 



2 2 



64tt 2 w„ A 4 64tt 2 i>„ A 4 



[1] I. Jack and D. R. T. Jones, Phys. Lett. B 457, 101 (1999) arXiv:hep-ph/9903365 . 

[2] I. Jack and D. R. T. Jones, Phys. Rev. D 61, 095002 (2000) arXiv:hep-ph/9909570 ; I. Jack, D. R. T. Jones and A. F. Kord, 

Phys. Lett. B 588 (2004) 127 [arXiv:hep-ph/0402045| . 
[3] M. S. Carena, S. Mrenna and C. E. M. Wagner, Phys. Rev. D 60 (1999) 075010 [arXiv:hep-ph/9808312l; M . S. Carena et 

al. [Higgs Working Group Collaboration] , |arXiv:hep-ph/0010338| U. Aglietti et al. , larXiv:hep- ph/0612i72] 
[4] M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, Nucl. Phys. B 453 (1995) 17 arXiv:hep-ph/9504378 ; M. S. Carena, 

S. Heinemeyer, C. E. M. Wagner and G. Weiglein, Eur. P hys. J. C 45 (2006) 79 7 |arXiv:hep-ph/0511023|; T. Hahn, 

S. Heinemeyer, F. Maltoni, G. Weiglein and S. Willenbrock, |arXiv:hep-ph/0607308| N. E. Adam et al, larXiv:0803.1154l 

[hep-ph]. 

[5] J. A. Casas, J. R. Espinosa, M. Quiros and A. Riotto, Nucl. Phys. B 436 (1995) 3 [Erratum-ibid. B 439 (1995) 
466] [arXiv:he p-ph/9407389 ; M. S. Carena, J. R. Espinosa, M. Quiros and C. E. M. Wagner, Phys. Lett. B 355, 
209 (1995) [arXiv:hep-ph/9504316| ; M. S. Carena, M. Quiros and C. E. M. Wagner, Nucl. Phys. B 461, 407 (1996) 

[arXiv:hep-plT/ 9508343 . 

[6] A. Sabanci, A. Hayreter and L. Solmaz, Phys. Lett. B 661 (2008) 154 [arXiv:0801. 20291 [hep-ph]]^ 

[7] M. B. Gavela, P. Hernandez, J. Orloff, O. Pene and C. Quimb ay, Nucl. Phys. B 430 (1994) 382 [arXiv:hep-ph/9406289] . 
[8] A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B 553 (1999) 3 [arXiv:h ep-ph/9902371 M. S. Carena, J. R. Ellis, S. Mrenna, 

A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B 659 (2003) 145 [arXiv:hep-ph/0211467| . 
[9] J. P. J. Hetherington, JHEP 0110, 024 (2001) [arXiv:hep-ph/0108206| ~ 
[10] H. E. Haber and J. D. Mason, Phys. Rev. D 77, 115011 (2008) [arXiv:0711.2890l [hep -ph]].; I. Antoniadis, E. Dudas, 

D. M. Ghilencea and P. Tziveloglou, AIP Conf. Proc. 1078, 175 (2009) |arXiv:0 809.4598 [hep-ph]]. 
[11] D . J. H. Chung, L. L Everett, G. L. Kane, S. F. King, J. D. Lykken and L. T. Wang, Phys. Rept. 407 (2005) 1 
[arXiv:hep-ph/0312378] . 

[12] L. Girardello and M. T. Grisaru, Nucl. Phys. B 194 (1982) 65j 

[13] J. Bagger an d E. Poppitz, Phy s. Rev. Lett. 71 (1993) 2380 [arXiv:hep-ph /93073 1 7| . 
[14] S. P. Martin, [arXiv:hep-ph/9709356| 

[15] H. E. Haber and R. Hempfling, Phys. Rev. Lett. 66 (1991) 1815; J. R. Espinosa and M. Quiros, Phys. Lett. B 266, 389 
(1991). 



[16] D. A. Demir, Phys. Rev. D 60 (1999) 055006 arXiv:hep-ph/9901389 ; A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B 

553 (1999) 3 arXiv:hep-ph/ 9902371 . 
[17] T . Ibrahim and P. Na th, Phys. Rev. D 63, 035009 (2001) [arXiv:hep-ph/0008237] , Phys. Rev. D 66 (2002) 015005 

[arXiv:hep-ph/0204092] . 
[18] D. A. Demir and L. L. Everett, Phys. Rev. D 69 (2004) 015008 [arXiv:hep-ph/0306240j . 

[19] R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH Collaboration and and], Phys. Lett. B 565, 

61 (2003) [arXiv:hep-ex/0306033] . 
[20] S. R. Coleman and E. J. Weinberg, Phys. Rev. D 7, 1888 (1973). 



10 



[21] M. A. Cakir, S. Mutlu and L. Solmaz, Phys. Rev. D 71, 115005 ( 2005) [arXiv:hep -ph/0501286]. 

[22] S. Y. Choi, M. Drees and J. S. Lee, Phys. Lett. B 481, 57 (2000) [arXiv: hep-ph/00022 87]. 

[23] S. W. Ham, S. H. Kim, S. K. OH and D. Son, Phys. Rev. D 76 (2007) 115013~ larXlv:0708.2755l [hep-ph]]. 

[24] A. De Roeck, V. A Khoze, A. D. Martin, R. Orava and M. G. Ryskin, Eur. Phys. J. C 25, 391 (2002) 
arXiv:hep-ph/0207042 ; M. S. Carena, S. Mrenna and C. E. M. Wagner, Phys. Rev. D 62, 055008 (2000) 
arXiv:hep-ph/9907422 ; M. S. Carena, H. E. Haber, H. E. Logan and S. Mrenna, Phys. Rev. D 65, 055005 (2002) 
[Erratum-ibid. D 65, 099902 (2002)] [arXiv:hep-ph/0106116| ;A. Duperrin, Eur. Phys. J. C 59, 297 (2009) |arXiv:0805.3624 
[hep-ex]]. 

[25] S. Abel, S. Khalil and O. Lebedev, Nucl. Phys. B 606, 151 (2001) arXiv:hep-ph/0103320 .; T. Ibrahim and P. Nath, Phys. 
Rev. D 58, 111301 (1998) [Erratum-ibid. D 60, 099902 (1999)] arXiv:hep-ph/9807501 . 



